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Abstract 

Form factors are derived for a model describing the coherent Joseph- 
son tunneling between two coupled Bose-Einstein condensates. This is 
achieved by studying the exact solution of the model in the framework 
of the algebraic Bethe ansatz. In this approach the form factors are 
expressed through determinant representations which are functions of 
the roots of the Bethe ansatz equations. 
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The experimental realisation of Bose-Einstein condensates in atomic alkali 
gases is generating a great deal of theoretical activity in order to understand 
the nature of coherent Josephson tunneling between coupled systems. A 
simple two-mode Hamiltonian which has been widely studied (see |1] for a 
review) takes the form 

H = \Kn 2 - - ^f(a\a 2 + a\ai). (1) 

o Z Z 

where a\,a 2 denote the single-particle boson creation operators, Ni = a\ai, 
N 2 = a 2 a 2 are number operators and n = Ni — N 2 is the relative particle 
number operator. Note that the total number operator N — Nx+N 2 provides 
a good quantum number since [H, N] = 0. 

While this model has been studied using a variety of techniques such as 
Gross- P it aevskii states Q, mean- field theory 0, quantum phase model || 
and an exact quantum phase model ||, it also admits an exact solution in 
the framework of the algebraic Bethe ansatz, given in H in the guise of the 
"discrete self-trapping dimer" model, which has been largely unexplored. In 
f|, we have used the exact solution to determine the asymptotic behaviour 
of the energy spectrum in the Fock (N << K/Sj) and Rabi (iV -1 >> K/Sj) 
regions. It was also shown that asymptotic thermodynamic properties can be 
deduced for the Fock region at low temperature and the Rabi region for all 
temperatures. Here, we will limit ourselves to the case Afj, = and continue 
our analysis of the Bethe ansatz solution to yield explicit exact form factors 
for the Josephson tunneling current 

J = i(a\a 2 — a 2 ai), (2) 

as well as the relative particle number n and n 2 , which are applicable for 
all couplings. This opens an avenue to explore the behaviour of the Hamil- 
tonian in the cross over between the Fock, Rabi and Josephson (iV -1 << 
K/Sj << N) regions. Recall that in Josephson's original work on macro- 
scopic superconductors [0], the tunneling current is a manifestation of the 
relative phase of the wave functions of two superconductors separated by an 
insulating barrier. As phase and particle number are canonically conjugate 
quantum variables, the quantum fluctuations of the relative particle num- 
ber are of primary importance in understanding the physics in the present 
model, given that there are technical difficulties which prevent a simple def- 
inition for the phase variable Jl], §. Our results for the form factors provide 
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an initial step towards the calculation of these fluctuations. The method 
that we will adopt follows that proposed for the form factors of the Bose gas 
with delta-function interactions and the one-dimensional Heisenberg model 
0, l0[ . This yields determinant representations for the form factors which 



are functions of the roots of the Bethe ansatz equations that arise from the 
exact solution. Furthermore, it is straightforward to express these results in 
a time-dependent form. 

First we will review the basic features of the exact solution of ([!]) via the 
algebraic Bethe ansatz, as discussed in ||. The theory of exactly solvable 
quantum systems in this setting relies on the existence of a solution of the 
Yang-Baxter equation 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v). (3) 

Here Rjk(u) denotes the matrix on V <8> V £g> V acting on the j-th and fc-th 
spaces and as the identity on the remaining space. The i?-matrix solution 
may be viewed as the structural constants for the Yang-Baxter algebra gen- 
erated by the monodromy matrix T(u) 



Ri2(u - u)Ti(«)T 2 (i;) = T 2 (v)T 1 (u)R 12 (u 



(4) 



For a given i?-matrix, there are a variety of realisations of the Yang-Baxter 
algebra. For the su{2) invariant i?-matrix, 



R(u) 



( 1 
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b{u) c{u) 

c{u) b{u) 








1 / 



(5) 



with the rational functions b{u) = u/(u + rj) and c(u) = f]/(u + rj), there is a 
realisation of the Yang-Baxter algebra in terms of canonical boson operators 



which reads 11 



L(u) 



u + r)N 



such that 



R 12 (u - v)L l (u)L 2 (v) = L 2 (w)Li(-u)i? 12 (M - v). 

The co-multiplication behind the Yang-Baxter algebra allows us to choose 
the monodromy matrix 

T{u) = Li{u + uj)L 2 {u — lo) 
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{u + u + r]Ni)(u- uj + 7]N 2 ) + <A a i (u + uj + r]N 1 )a 2 + r] 1 a 1 
(u — u + r)N 2 )a\ + rf x a} 2 a\a 2 + rf 2 

(6) 

Defining the transfer matrix through t(u) = tr (T(u)), it follows from ([|) that 
the transfer matrices commute for different values of the spectral parameter; 
viz. 

[t(u),t(v)] = Vu,v. 
In the present case, we have explicitly 

t{u) = u 2 + ut]N + rfNiN 2 + rjuh + a) 2 a>\ + a\a 2 + rj~ 2 — uj 2 . 

Then 

dt. 

T]N 



m = * 

du 



u=0 

and it is easy to verify that the Hamiltonian is related with the transfer 
matrix t(u) by 

H=-k (t(u) - ^(t'(0)) 2 - ut'(0) - if 2 + uj 2 - u 2 ^ , 
where the following identification has been made for the coupling constants 

K K7] 2 

T ~ ~2~' 

— = -KTJUJ, 

2 

The solution of (|T|) via the algebraic Bethe ansatz is obtained by utilizing 
the commutation relations of the Yang-Baxter algebra. Setting 

we have from the defining relations (§) that 

[A{u),A(v)] = [D(u),D(v)] = 0, 
[B(u),B(v)} = [C(u),C(v)} = 0, 

A(u)C(v) = U ~ V + V C(v)A(u) ^— C(u)A(v), 

u — v u — V 

D(u)C(v) = U ~ V ~ V C(v)D(u) + -H—C(u)D(v). (8) 
u — v u — V 
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An explicit representation of (|7]) is obtained from (^) with the identification 

A(u) = (u + lu + rjNj (u-uj + r]N 2 ) + a\ai 

B(u) = (u + uj + r]Ni)a 2 + rf l a\ 

C(u) = (u-uj + r)N 2 )a\ + rf l a\ 

D{u) = a\a 2 + rf 2 . 



Choosing the state |0) = |0) 1 |0) 2 as the pseudovacuum, the eigenvalues a(u) 
and d{u) of A{u) and D{u) on |0) are 

a{u) = {u + oj){u — lj), 
d(u) = r^ 2 . 

Next choose the Bethe state 

AT 



\V 



) = \v 1 ,...,v N ) = Y[C(v i )\0). (9) 



i=l 



Note that because [C(u), C(v)) = 0, the ordering is not important in (Jfty. It 
is also clear that these states are eigenstates of N with eigenvalue N. The 
approach of the algebraic Bethe ansatz is to use the relations (|8]) to determine 
the action of t(u) on \v). The result is 

t(u) \v) = A(u, v) \v) 

V a ( v i) TT v i ~ v j + V 1 i \ 

2^ 11 \ v i, -Vi-i,u, v i+1 , ...,v N ) 

fu-Vif^ Vi-Vj J 
\ a U ~ Vi f£ Vi - Vj J 

where 

\( -A f \T\ U ~ Vi + 7] i M \TI U ~ Vi ~ 1] (^\ 
Am, u) = a(u) h . (11) 

l=1 u-Vi i=1 U-Vi 

The above shows that \v) becomes an eigenstate of the transfer matrix with 
eigenvalue ( |TT| ) whenever the Bethe ansatz equations (BAE) 

N 

^ Vi -Vj+T] 
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are satisfied. As N is the total number of bosons, we expect N + 1 solutions 
of the BAE. Note that in the derivation of the BAE it is required that 
Vi 7^ Vj Vi, j. For example, the solution 



of (0) is invalid, except when N = 1. This is a result of the Pauli Principle 



for Bethe ansatz solvable models as developed in |L2] for the Bose gas. We 
will not reproduce the proof for the present CclSC, clS it follows essentially the 
same lines as ||12||. For a given valid solution of the BAE, the energy of the 



Hamiltonian is obtained from the transfer matrix eigenvalues ( |TT| ) and reads 

( N n 2 N 2 \ 
E(v) = — k I rj^Wrfivj - u + r])(vj + u) - tjujN - r)~ 2 J . (13) 

A consequence of the above construction is that form factors, and conse- 
quently correlation functions, can be computed for this model. The method 
we employ follows that used for the Bose gas and one-dimenisonal Heisenberg 
chain || [I0[] . The reason we can adopt this approach is because the solution 
of these models and ([I]) are based on the same i?-matrix (||). Below we will 
restrict our attention to uj = (for reasons which will become clear) and give 
explicit results for n, n 2 and J . 

A curiosity of this model is that the representation of the Yang-Baxter 
algebra is non-unitary; viz. 

C\u) ^B{u). 

In the case when uj = it is, however, equivalent to a unitary representation 
since in this instance we have 

C\u) = P.B(u).P 

where P is the permutation operator. The permutation operator is defined 
by the action on the Fock basis 

P(aly(al) k \0) = (a\) k (aty\0). 

Note that for uj = (which will be assumed hereafter), we have 

[P, H] = 
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which means that the energy eigenstates are also eigenstates of P. Moreover, 
P 2 = I shows that P has eigenvalues ±1. 

There exists a formula due to Slavnov [13| for the scalar product of states 
obtained via the algebraic Bethe ansatz for the i?-matrix (H), which when 
applied to this model is 

S(v:u) = (0\B(v 1 )...B(v N )C(u 1 )...C(u N )\0) 

det F(u : v) 

where 

v -l I N N 

Kb = — \ u 2 b Y[ (v m - u b - rj) - J] 2 Yl ( V m -Ub + rj) 



u b - v a 



the parameters {v a } satisfy the Bethe ansatz equations fljjD and {uj} are 
arbitrary. Note that when \u) = \v) we need to take a limit for the diagonal 
entries which gives 



N N N 



F aa = -2r/ 2 v a Y[(vi -v a -r])-r) 4 Y[(vi - v a + rj) 

i i 



^ a V m -V a + T] 



N N -2 2 

+n^-« B -i?)E _ _ • 

I m^a Vm Va V 

Consider 

S(v:v) = (0\B(vi)...B(v N )C(v 1 )...C(v N )\0) 

= (0\PC\v l )...C^v N )PC(v l )...C(v N )\0 

= e(v) (v\v) 

where e(v) = ±1 is the eigenvalue of P; viz. 

P \v) = e(v) \v) . 

This quantity can be determined to be given by 

N 

e{y) = l[r)Vj 
j 



Hence, from the Slavnov formula, the norms of the eigenstates 

\\v\\ = \v\v) 



\S(v:v)\ 1/2 



are obtained directly. 
Let us define 



X = A(0) - D(0) = r?N x N 2 rf 2 
where J is defined by (0). In analogy with ([LT]) we find 

X\u) 

= 0(u) \u) 

* r]a{ Ui ) ILui-Uj + r] 
+ 2^ 11 - — Fl' ■■■ u i-li °) -,Un) 

" rjd(ui) JLui-uj-T) , 
+ 2^ 11 _ \ u u —Ui-i, 0, Ui+i, u N ) 



where 



N AT i 

i=l u * i=l u i 



= -v 11 

t=i u * 

Using the Slavnov formula we can calculate the form factors for x- Suppose 
that \v) is an eigenstate of the Hamiltonian. Then 

(v\x\u) 
= 6{u)e{y)S{v : u) 

. A va(ui) SL ui - Uj + 7] 

+ 2^ 11 — e(v)S(v : ui, ...Mi_i,0,M i+ i, ...,u N ) 

+ 2^ 11 _ e{v)S{v : Mi, ...Ui-i,Q,Ui+i, ...,u N ). 
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Now we let \u) be an eigenstate of the Hamiltonian (J!]). The above formula 
can be simplified considerably (cf. fT0|l ) 



(v\x\u) 



~r] N ~ 2 Ua(Va + V) 



det (F(0: u) - 2r]' 3 e(u)e(v)Q(v : u) 




Uj>k(u k - Uj) Y[ a< p{vp - v a ) 




This is our main result. We remark that because the basis states are also 
Hamiltonian eigenstates, it is straightforward to write down the time-dependent 
form factors 



where the energies E(v) are given by (|T3|). Remarkably, from eqs. ([14], [15]) 
all the time-dependent form factors for n, h 2 and J can be obtained. This 
is achieved by exploiting the symmetry of the Hamiltonian under P, which 
we will now explain. 

The following corollary is easily proved. If e(v) ^ e(u) then 



(v\x(t)\u) = exp(-it (E(v) - E(u))) (v\x\u) 



(15) 




If e(v) = e{u) then 

(v\J\u) = 0. 



The result follows from the observation 



PiViiV 2 
PJ 



N!N 2 P, 
-JP 



We now find that 



(yiNx^lu) = "if 2 (v\x(t)\u) + rf A (v\u) 



if e(v) = e(u), and is zero otherwise. Also 



(v\J\u) 



i (v\x(t)\u) 



9 



if e(v) 7^ e(u), and is zero otherwise. 

The above shows that the form factors for J are obtained directly from 
those of x(t). Those for n 2 also follow, since we have n 2 = N 2 — AN1N2 and 
the states (P) are automatically eigenstates of N with eigenvalue N. Thus 



(v\n 2 \u) = N 2 (v\u) - 4 (y\N x N 2 \u) 



To obtain the form factors for n, we use the fact that J is the time derivative 
of n, so 



which gives 



J=-[n, H] 

K 



(v\n\u) = m\u) • 



In principle, the correlation functions 

mom 



(*l*> 



where 9 = n, h 2 or J ', and is an arbitrary state, can be expressed in 
terms of the form factors through completeness relations. In particular, for 
a given the quantum fluctuations of the relative number operator 



can be computed from these results. 

In conclusion, we have shown that the algebraic Bethe ansatz solution 
of ([]]) provides a means to calculate form factors, and in turn correlation 
functions, for the model. We have yielded explicit results for the case Afi = 0. 
An outstanding task is to extend these results to the case A/x 7^ 0, which 
presents a challenging problem because of the non-unitary representation 
of the Yang-Baxter algebra in this instance, and the breaking of symmetry 
under P. 
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